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Abstract 

We consider the linear elliptic eqnation —div(oVu) = / on some bounded domain D, where 
a has the affine form a = a{y) = a + J2j>i Vji’j some parameter vector y = (yj)j>i G U = 
[— 1 , 1 ]^. We study the summability properties of polynomial expansions of the solution map 
y I—)■ u{y) G V := Hq{D). We consider both Taylor series and Legendre series. Previous results 
|S] show that, under a uniform ellipticity assuption, for any 0 < p < 1, the C summability of the 
(ll'0ilU°°)j>i implies the C summability of the C-norms of the Taylor or Legendre coefficients. 
Such results ensure convergence rates of polynomial approximations obtained by best n- 
term truncation of such series, with s = ^ — 1 in L°°{U, C) or s = ^ ^ in Lp‘{U, V). In this 

paper we considerably improve these results by providing sufficient conditions of summability 
of the coefficient C-norm sequences expressed in terms of the pointwise summability properties 
of the {\fpj\)j>i- The approach in the present paper strongly differs from that of [S], which is 
based on individual estimates of the coefficient norms obtained by the Cauchy formula applied 
to a holomorphic extension of the solution map. Here, we use weighted summability estimates, 
obtained by real-variable arguments. While the obtained results imply those of [S] as a particular 
case, they lead to a refined analysis which takes into account the amount of overlap between the 
supports of the ijjj. For instance, in the case of disjoint supports, these results imply that for 
all 0 < p < 2, the C summability of the coefficient Wnorm sequences follows from the weaker 
assumption that (||■^/'ilU“)i>l is summable for q = q(j}) := > p. We provide a simple 

analytic example showing that this result is in general optimal and illustrate our findings by 
numerical experiments. The analysis in the present paper applies to other types of linear PDEs 
with similar affine parametrization of the coefficients, and to more general Jacobi polynomial 
expansions. 
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1 Introduction 


1.1 Elliptic PDEs with afRne parameter dependence 

The numerical treatment of parametric and stochastic partial differential equations was initiated in 
the 1990s, see |10l [TTl [T3l [T6] for general references. It has recently drawn much attention in the 
case where the number of involved parameters is very large Ei, or countably infinite [I10I1]. 

In this paper, we are interested in the approximation of the elliptic parametric PDE of the form 


— div(aVn) = /, 


( 1 ) 


set on a bounded Lipschitz domain H C (where in typical applications, d = 1,2,3), with 
homogeneous Dirichlet boundary conditions. We consider coefficients having the affine form 


a = a{y) = a + 

i>i 


( 2 ) 


where 

ye [/ = [-!, if. (3) 

The functions a and ipj are assumed to be in L°°{D). 

Well-posedness of this problem in V = Hq{D) is ensured for all y G U by the so-called uniform 
ellipticity assumption, further referred to as (UEA), 


<d-r. (UEA) 

i>i 

for some r > 0, which is equivalent to r < a{y) < R for all y G U for some 0 < r < i? < oo. This 
assumption ensures the uniform boundedness of the solution map y i—)> u{y) from U to V. 

We are interested in polynomial approximations of the solution map. Such approximations have 
the form 

Uk{y) ■=^Vuy'', (4) 

where A C T" is a hnite set of (multi-)indices ly = {iyj)j>i G T and y^ = nj>i Vj- the inhnite 
dimensional setting, the index set T denotes the (countable) set of all sequences of nonnegative 
integers which are finitely supported (i.e. those sequence for which only hnitely many terms are 
nonzero). Note that the polynomial coefficients Vi, are functions in V. Such approximations can be 
obtained by truncation of inhnite polynomial expansions. 

Two relevant such expansions are: 


1. The Taylor series 

tt = ^ Uy’', U := ^d''u{0), 

vGT 


(5) 
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where we use the customary multi-index notations 


y'' = Y[ Vj' and 

i>i i>i 

with the convention that 0! = 1. 


( 6 ) 


2. The Legendre series 

u = ^ ^ Uj/Liy, Uy = Ljyijj) = Luj iyj ), (7) 

v&T j>l 

where the univariate Legendre polynomials are normalized in L^([—1,1], y), and (•, •) 
denotes the inner product in 1/,/r), where /i = (^j>i ^ is the uniform probability 

measure. 


Given such expansions, one natural way of constructing a polynomial approximation is by best 
n-term truncation, that is, setting 

ul{y) ■■= Uy’", (8) 

or 

Uniy) ■= Y ( 9 ) 

where and are the index sets corresponding to the n largest ||tt/||v" or ||Ui/||\/, respectively. 
The convergence rates of such approximations are governed by the P’ summability properties of 
the sequences for p < 1 and (||n,y||v")i/eJ' for p <2. Indeed, by standard application of 

Stechkin’s lemma [9], such summability properties imply error estimates of the form 


and 


W ~ ^ C{n -I- 1) 


“ Un\\L'^{uy,iJL) < C{n -h !)■ 


s - 1, 
P 


C := 




y\\V )v(iT\\£P^ 


\Uv\\v)u&t\ 


( 10 ) 


( 11 ) 


see [8]. Let us mention that the case of hnitely many variables, that is, when -i/’j = 0 for j > d, 
typically leads to exponential convergence rates of the form exp{—cn^^^). Such rates have been 
obtained in mm by application of the above Stechkin estimates for all 0 < p < 1 and tuning the 
value of p with that of n. They can also be directly obtained through the available upper bounds 
on the coefficients ||ti/||v" or ||ui/||\/, as shown in [15], which brings certain improvements over the 
previous approach. In the present paper we focus our attention on the inhnite-dimensional case. 

The above approximation estimates have general implications on the potential performance of 
other model reduction techniques. Indeed we observe that u^{y) and u^{y) belong to hxed n- 
dimensional subspaces of V, spanned by the {tu)y^jsT or It thus follows from (10) that 

the best re-dimensional model reduction error in the uniform sense, given by the Kolmogorov re-width 


dni-M.) := inf sup min ||u — rc||y, M := u{U) = {u{y) : y£U}, (12) 

dim(V,i)=n ^g_A4 w&Vn 
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decays at least as fast as n Likewise, it follows from ( [TT| that if y is uniformly distributed in U, 
the best n-dimensional model reduction error in the mean-square sense 

inf E{mm \\u{y)-w\\l), (13) 

dim{l4,)=n wGVn 

attained when Vn is the span of the n first ^-principal components, decays at least as fast as 
^_ 2 r above estimates govern the convergence rate of reduced basis (RB) methods and proper 
orthogonal decomposition (POD) methods, respectively, see |7]- 


1.2 Existing results 

Summability results for the sequences (||ti,||y)jyg_ 7 r and (||M,y||y)i,g_ 7 r have been established in | 8 ] 
under the so-called uniform ellipticity assumption 

'^\ipj\<a-r. (UEA) 

i>i 

for some r > 0, which is equivalent to r < a{y) < R for all y € U for some 0 < r < i? < oo. This 
assumption ensures the uniform boundedness of the solution map y i—>■ u{y) from U to V. These 
results can be summarized as follows. 


Theorem 1.1 Assume that (UEA) holds. For any 0 < p < 1, if (||V'j||L°°)j>i £ .^^(N) then 
{\\U\\v)ugT e and {\K\\v)u&t G 

These results have been extended to a large range of linear or nonlinear parametric PDEs mm 
where y is again ranging in the infinite dimensional box U. They strongly rely on the holomorphy 
of the solution map y i—>■ u{y) in each variable yj. While they yield provable algebraic approxima¬ 
tion rates for parametric PDEs in the infinite dimensional framework, these results are not fully 
satisfactory for several reasons. 

First, they are confined to the case 0 < p < 1. Therefore, no summability property of the Taylor 
or Legendre coefficients can be deduced from these results in cases where (||V’ilU°°)i>i does not 
have this type of summability. Note also that, in the case of Legendre series, convergence rates 
of best re-term approximation can be derived as soon as (||rei/||\/)i/gj- G for some p < 2, in 


view of (11). A legitimate objective is therefore to identify conditions on {il>j)j>i that govern the 
iP summability of (||rei,||\/)j/eJ' the regime 1 < p < 2 . 

Second and more important, the summability conditions imposed on the ||^j||ioo in this result 
becomes quite strong and artificial in the case where the supports of these functions do not overlap 
too much. As a relevant example, consider the case where the are a wavelet basis on the 

domain D. In this case it is more natural to denote the elements of such bases by fjx, where A 
concatenates the space and scale indices, following the usual terminology such as in |^, and using 
the the notation / = |A| for the scale level. At a given scale level I > 0, there are 0(2'^^) wavelets 
and their supports have finite overlap in the sense that, for all x € D, 


if{X : |A| = Z and if\{x) 7 ^^ 0} < M, 


(14) 
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for some fixed M independent of I, where denotes the cardinality of the set S. It is well known 
that the geometric rate of decay of the wavelet contributions as the scale level grows reflects the 
amount of smoothness in the expansion. It is thus natural to study the situation where 

= Cl := C2-^\ |A|=/, (15) 

for some given a > 0. In the affine case, assuming that a and the tpx are sufficiently smooth, this 
means that the diffusion coefficient a is uniformly bounded in the Holder space C^{D) independently 
of y. In this case, the uniform ellipticity assumption is ensured provided that 

M Cl < ajnin — f, where amin := ess inf/j a > 0, (16) 

l>0 

which holds if C < M~^{aram ~ ^)(1 ~ 2““) for the above constant C. Note that we can take a > 0 
arbitrarily small, up to taking C small enough. The uniform ellipticity assumption ensures that 
the solution map belongs to L°°{U,V) and therefore to L‘^{U,V, y), and thus that the Legendre 
coefficient sequence (||rii/||v")j/eJ' belongs to However, if we want to use the above mentioned 

results to prove summability of this sequence for smaller values of p, we are enforced to assume 
that the are summable over all indices, which equivalently means that 

< oo, (17) 

l>0 

and therefore that a > d. This is highly unsatisfactory since it excludes diffusion coefficients with 
low order of smoothness, especially when the spatial dimension d is large. 


1.3 Main results and outline 


The above example reveals a gap in the currently available analysis: summability of (||?;,y||y)[/g^ 

can be obtained under mild assumptions, while proving summability for p < 2 by the existing 
results immediately imposes much stronger assumptions (in the sense of the required decay of 
Il'i/ijllioo as j —?• oo). In this paper, we propose a new analysis which allows us to remove this gap. 
The main results of the present paper can be summarized as follows. 


Theorem 1.2 Assume that ( UEA) holds and that for 0 < p < 2 and q = q{p) := there exists 
a sequence p = {pj)j>i with pj > 1, such that 


(18) 

i>i 


for some s > 0, and sueh that iyp- belongs to .^'^(N). Then the sequences cmd 

{\WAv) belong to 


This theorem constitutes a particularly strong improvement over Theorem |l.l| when the supports 
of the functions {ifj)j>i do not overlap much. For example, in the case where these supports are 
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disjoint we can take weights such that p- ^ ~ ||^j||x,oo as ||V’jHL<=o tends to 0. Therefore, in this case, 
for all 0 < p < 2 the (P summability of and (||rijy||y),ygj- follows from the assumption 


that 


'j\\L°°)j>i is summable for q = q{p) := 2 ^- Note that q{p) > p and that 

lim q(p) = + 00 , 

p—>-2 


(19) 


which shows that almost no decay of {\\'tpj\\L°°)j>i is required as p gets closer to 2. Similar improve¬ 
ments can be obtained for other types of families such as wavelets. 

Let us also mention that, while we focus on the diffusion equation, inspection of proofs reveals 
that the main results can be extended to other types of linear elliptic or parabolic PDEs with similar 
affine dependence of the coefficients. Specific examples are given at the end of §2. 


The remainder of this paper is organized as follows. In §2, we establish Theorem 1.2 in the case 
of Taylor coefficients. The approach strongly differs from that used in for proving Theorem o 
which is based on establishing individual estimates of the ||ti/||y by Cauchy’s formula applied to 
to a holomorphic extension of the solution map. Instead, we use weighted summability estimates, 
obtained by real-variable arguments. It is worth mentioning that these estimates also imply the 
individual estimates from j8]. 

In §3, we establish Theorem 1.2 in the case of Legendre coefficients. The approach is again 


based on weighted summability estimates, obtained by combining the ideas developed for the Taylor 
coefficients with Rodrigues’ formula. Furthermore, we give an extension of our results to Jacobi 
polynomial expansions by an analogous argument. 

We give in §4 several examples of applications, corresponding to different types of support 
properties for the {ipj)j>i and we discuss in each case the improvements over Theorem In 

particular, we show that our results are sharp in the case of disjoint or finitely overlapping supports, 
in the sense that summability of (HV’jwith q := q{p) is generally necessary to achieve 
iP summability of (||t[/||y)i/gj-. While we focus on the infinite-dimensional framework, the obtained 
results can also be used in the case of finitely many variables. In particular, one may use them to 
obtain improved exponential convergence rates, by following either a Stechkin-based approach as in 
mm or a more direct approach as in |15| . 

Finally, we give in §5 the results of several numerical tests evaluating the sharpness of the 
obtained results for the various types of (V’j)j>i considered in §4. 


2 Summability of Taylor coefficients 

We start from an analysis of Taylor coefficients, where we use the following alternative expression 
of the uniform ellipticity assumption, further referred to as (UFA*), 


E7>1 IV’il 


< 1 . 


(UFA’' 


L°° 
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Indeed, if (UEA) holds with some r > 0, then we also have 


Sj>i IV’il 


< 


a — r 


< 1 - 




< 1, 


( 20 ) 


and it is also easily checked that (UEA*) implies (UEA) for a certain r. 


Similar argnments as in [H] show that nnder snch an assnmption, the partial derivatives d'^u{y) 
are well defined for each v € as elements of V for each y G U. They can be compnted by applying 
the operator in the y variable to the variational formnlation 


Id 


a{y)Vu{y)-Vvdx = {f,v)v',v, v G V, 


( 21 ) 


which gives, for all ^ 0, 


ID 


a{y)Vu{y) ■ VV dx = — ^ Uj j u{y) ■ Vv dx, v G V, 

jGsupp u 


ID 


where 


Bj — (0, . . . , 0, 1, 0, . . . ) — (5jj)i>l 

is the Kronecker 5 seqnence with 1 at position j. For the Taylor coefficients, this yields 


aVt,y • S/v dx = — 


ID 


/ il)jVty-ei • Vn dx, V gV, 


jEsupp u 


when u ^ Q and 


/ aVto ■ Vv dx = {f,v)v',v, vGV. 
Jd 

We shall make nse of the norm defined by 

||n||a := / a|Vr;|^ dx. 

Jd 

This norm is obvionsly eqnivalent to the U-norm with 

dminiblly < ll^^lll < ||a||Loo||n||y. 


( 22 ) 


(23) 


(24) 


(25) 


(26) 


(27) 


The following results shows that under (UEA*), the energy of Taylor coefficients decays geometri¬ 
cally with the total order of differentiation. Its proofs uses a technique introduced in |5]. Here we 
use the standard notation 

(28) 


i>i 


to denote the total order. 
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Lemma 2.1 If (UEA*) holds, then for a := < 1 with 


6 = 


Z^j>i IV'il 


< 1, 


we have for all k > 1, 


\i>\=k—l 


.1,1 


Proof: For u £ we define 

d,y := / d\Vt,^\‘^dx, := / \'il)j\\Vt,j\^ dx . 

Jd ’ Jd 

Then by ( |UEA*[ ), 

^ ^ di/^j ^ Odi/. 
i>i 

Fnrthermore, since for i' ^ 0 the Taylor coefficients satisfy 


aVt,y • Vn dx = — 


ID 


y~] / fjjVtu-ej ■ Vn dx, V £V, 


jEsupp u 

we also have 

d.< Y. myu-e,\£tAdx<\ Y {du-ej ,j + ) ) 

jSsuppi^ ^ jSsupp;/ 

where we have nsed Yonng’s ineqnality. Thns, by 

e 


\ E 


jEsupp V 


Snmming over \v\ = k, we obtain 


1-? 

2 


Y1 du-e,,j = l Y1 

\u\=k \u\=k jeswppu \u\=k-l j>l 


and nsing again (32), we arrive at (30). 


(29) 

(30) 

(31) 

(32) 

(33) 

(34) 

(35) 

(36) 

□ 


Note that as an immediate conseqnence of (30) 

Ewiis 


2 - 


ueD 


I-a 2-26 


a<00, 


and therefore by (|27[) and the Lax-Milgram stability estimate for tp = n(0), 


Ell* 

veT 


i/||y E 


(2 — 0) IIoIIloo 
(2 - 20)0111111 


V ^ 


(2 — 0) ||a||£,oo 

(2-20)a^iii 


V 


/ < oo. 


(37) 


(38) 


We are now ready to state the main resnlt of this section which is a direct conseqnence of the above 
observations. 
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Theorem 2.1 Assume that for some sequenee p = {pj)j>i of positive weights, we have the weighted 
uniform ellipticity assumption 

Ej>i Pj\^j\ 
d 

Then the sequence {p'^\\tv\\v)u&F is summable, with 

Y,iP''\Mvf<C<oo, (40) 

v^T 


< 1 . 


(39) 


where 


C = C{d,f,6) 


(2 — (5) ||a||ioo 

(2 - 



(41) 


Proof: We observe that the weighted UEA is equivalent to (UEA*) for the rescaled coefficient 
Opiu) ■= o,{Dpy), where DpU := {pjyj)j>i- So we obtain with Lemma 2.1 and (38) that 



'y ^ II ^Pit'll ^ 

(42) 




where 

tp,v = \d''up{tf) = p'^u, Up{y) = u{Dpy). 

(43) 


I/! 


The result follows. 


□ 


As a consequence we obtain the following summability result. 


Corollary 2.1 If for some sequence p = {pj)j>i with pj > 1, j € N, we have the weighted uniform 
ellipticity assumption (39) and if the sequence belongs to .^'^(N) with q = q{p) := for 

some p <2, then the sequence (||d| y) is P’ summable. 


Proof: By Holder’s inequality, 

Moreover 


2 

i2||y 


i/eJ' 




p/2 


(E- 


X (2-p)/2 


2-p- 


-'=n(E^-, 

u&T j>l k=0 


J ~ 


i>i 


(44) 


(45) 


where the latter product converges precisely when {p^ )j>i E The statement thus follows from 


Theorem 12.11 


□ 


Remark 2.1 As a trivial consequence of Theorem 2.1 using the fact that the norm dominates 
the norm, we also retrieve the estimate 


tv\\v < Cp ^, 


(46) 
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for any sequence p such that (39) holds, where C is the square root of the constant in (41 )■ This 


estimate was established in m, with a different constant, by a complex variable argument, namely 
invoking the holomorphy of y u{y) on polydiscs of the form ®j>i{\zj\ < pj}. One advantage 
of this individual estimate is that one may choose to optimize over all possible sequences p, which 
yields 

\\tu\\v < C inf p~’', ( 47 ) 

where the infimum is taken over all sequences p = {pj)j>i of numbers larger than 1 such that 
(39) holds. The proof of Theorem \n\ in m for the Taylor coefficients is based on using the above 
estimate, which amounts to selecting a different sequence p = p{v) for each v d T. However, we 
show in §4 that in several relevant cases better results can be obtained by using Corollary\2.1\ 


The above analysis exploits the particular structure of the linear diffusion equation and affine 


dependence of a in the variables yj, in particular when deriving the recursion equation (24) which 
subsequently leads to the estimate of Lemma 2.1 It can be extended in a natural manner to 


other types of equations with similar properties. Here are two simple examples. 

1. The fourth order equation A (a Art) = / set on a Lipschitz domain D with homogeneous 


boundary condition tt = = 0 on dD and a having the affine form 

relevant space is V := Hq{D), the closure of T){D) in and the recursion is 


In this case the 


aAt,,Av dx = — 


ID 


E 

jGsupp 12 


fjAty-e.Avdx, V G V. 


( 48 ) 


ID 


By taking v = ty we obtain Lemma 2.1 for the norm ||w||a := (Jp a\Avf')^^‘^, which is equiv¬ 
alent to ||f||y = ||At;||j;^ 2 . 

2. The parabolic equation dfU — div(aVtt) = / set on ]0,T[xD with homogeneous boundary 
condition u = 0 on ]0,T[xdD and initial condition u^t=o = uo ^ this case, the 

relevant solution space is V = L^(]0, T[; L7 q (H)) n H^(]0,T[; H~^(D)). We refer to [H] for 
the corresponding space-time variational formulation and the analysis of its well-posedness. 
For y 7 ^ 0 we find that ty satisfies a parabolic PDF with initial value tyj^^Q = 0 and recursive 
variational form given by 

f f {dttyV+ aS/ty ■Vv)dxdt = — f f fjVty-ej-Vvdx, v G V. ( 49 ) 

Jo Jd Jd 


^Esupp u 


By taking v = ty we obtain Lemma 2.1 for the norm HwHa := (Jq Jd alVup)^'^^. Note that 


this norm is not equivalent to the full norm of V but only to that oiV = L^(]0, T[; Hq{D)). 

3. Similar results hold when a is of symmetric tensor type, with d non-degenerate uniformly over 

1 /Q 

D. We obtain Lemma 2.1 for the norm ||f||a := (/£)h|Vr;p) , under a uniform ellipticity 

assumption taking the following form, 

^i>i 


^ill2 


Amin(fl) 


< 1 . 


(50) 
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analogous to with HV’jlb and Amin(o) denoting the functions x i—)• ||V'j(x )||2 (with || • ||2 

the spectral norm) and x Amin(h). Note that this assumption is equivalent to p9[ ), that is, 
to (UEA*) or (UEA) in the scalar case. Likewise we obtain the subsequent results under the 


assumption that the sequence {pj)j>i satishes 


E,>iPjllV’ill2 


'^min (®) 




< 1 . 


3 Summability of Legendre coefficients 

in this section, we show that the summability properties of Corollary |2 .1 1 hold also for the Legendre 
coefficients of u. 

Theorem 3.1 If for some sequence p = {pj)j>i with pj > 1, j £ N, we have the weighted uniform 
ellipticity assumption 

Ylj>i Pj\'4’j\ 


S := 


< 1 , 


(51) 


then with := nj>i sequence p'^\\Uy\\v)v&T is ^ summable, that is, 

'^{a~^p''\\uv\\v)‘^ <C <oo, (52) 




where 


C = C{a,f,6) := 


(2-5)il + 6)\\-a\\U\frv, 


(53) 


2(1 - 

Proof: For y,z £U, we set TyZ := [yj + (1 — \yj\)pjZj )Then for Wy{z) := u{Tyz), we have 

irw,( 0 ) = (n(l - ls,l)''')p"8‘'o(!/). (54) 

i>i 


Let us hx y G C and set Oy := a{y) = a + X]j>i 'Uj'^Pj and ij^yj := (1 — \yj\)pj'4’j- Then Wy is the 
solution of 


— div 




(55) 


i>i 


Applying Lemma 2.1 with the modihed y-dependent coefficients in (55), for the Taylor coefficients 


ty,u '■= (z^O ^d'^Wy{0) of Wy we obtain 


Py.i'llay ^ Pi/.i'llayi ZJ. 

\u\=k L|=fc—1 




^~2-ey’ y 


X]i>i \'^y,j\ 


(56) 


Since Pj > I, 


dy < 



< 

Ej>i Pj h/’i 


L°° 

a 


= (5 < 1 , 


(57) 
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and thus ay < S/{2 — S) < 1. Consequently, as in (38), 


— W^y ®IU°° Py,i'llay — W^y 


v^T 


v&F 


-II ^ ^ 11+ l|2 


(58) 


where 


We also have 


||t?;,o|lay < ||ay||L°°||«y Wl^WJWV 


(59) 


ll^yII l*^ ll^y^ ||^||l°° 


i>i 


a — 


i>i 


-3 


S (60) 

L- (1 - 5)3a^i^ 


by (51). Altogether, we obtain 


V lit lP<a“^Vlh ||2 ^ (2-<^)(l + (^)l|Q|lioo||/lly/ 

/ , \ry,i'\\v ^ ^min / ^ Iry.t^lla — — O < oo. 


i+eJ' 




{2-26){l-6r-ai^ 


(61) 


+ 1^2 




(y -i)'^^ ■ 


= / u{y)L^{y) d^{y) = / 'a(y) TT 5; 
Ju Ju t-C, 


j>i 




(y2 - l)-o dy{y) 


n \/2)^) I ^9"n(y) [] 

.,'^1 JU ' 


J>1 


i>i 


-^,— dy{y). 


(62) 


With the present normalization as in ([^, the Legendre polynomials satisfy the Rodrigues’ for¬ 
mula 

^.(y) = n^% 

i>i 

As a consequence, for the Legendre coefficients of u we obtain 

V^2zzj + 1 


(63) 


Hence, by (54) and (58), 




u&r j>l 


v^T 


W 


fn{y) ^ n 


2 -TT (1 - v]?''" 


i>i 


221- 


dp(y) 


Y.f’-'' ) 7 «''“wLn(i-iK 


/rr --- I U\ 


i>i 




- / ^dfiiy)= [ ^\\tyAvdy{y) <C, 


which completes the proof. 


□ 


Corollary 3.1 If for some sequence p = {pj)j>i with pj > 1, j G N, we have the weighted uniform 
ellipticity assumption (Hg, if the sequence (p^ belongs to with q = for a p < 2, then 

the sequence (lliti+lly)^,^^- is summable. 
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Proof: We obtain the statement from Theorem |3.1| by proceeding exactly as in the proof of Corollary 
12.11 In this case we need 

(65) 


■i>i 


which, since pj > 1, holds precisely when {p- ^)j>i G since 


E+1)"'"=n(E + 

j>l j>l k=0 


( 66 ) 


and since, by nsing the fact that \\{pj )j>i||£oo < 1, we find that the snm in each factor of the above 
prodnct converges and is bonnded by 1 + Cp~^ for some fixed C. □ 


Remark 3.1 Similar to the case of Taylor coefficients, we can also derive from Theorem 3.1 the 
individual estimate 

\Wv\\v < Ca„p~'^, (67) 

which is very similar to, yet slightly better than, the one established in m by complex variable 
arguments. 


Remark 3.2 As a consequence of Stechkin’s lemma, under the assumptions of Corollary \3.1\ we 
find that the best n-term approximation polynomials 


u’d := 


•- ^ ^ UyLy, 

obtained by retaining the indices of the n largest ||ttjy||y; satisfy the estimate 


( 68 ) 


< n ^ (69) 

where ^ ^ ~ 5 = ^- There is, however, a more direct and constructive way of retrieving this 

convergence rate, namely taking instead to be the set of indices corresponding to the n smallest 
values of the weights Wy := a~^p'^ which appear in |5^. We then directly obtain that 

\i/2 


\U - Un\\L2{U,V,u) ^ sup W,, 

U^An 




< 3* 

~ “n+l! 


(70) 




where (d)i)n>i is the decreasing rearrangement of the sequence (rCj^ As seen in the proof of 

Corollary 3.1. this sequence belongs to C^{T), which implies that d* < n~^ with r := ^. 

Remark 3.3 The results of this section can be generalized to other families of orthogonal polyno¬ 
mials satisfying a suitable Rodrigues’ formula. For instance, when the uniform measure on U is 
replaced by a tensor product beta measure 


dp{y) = 


TlyOLj + /3j + 2 ) 




^ 2«.+/3i+ir(a^- + + 1) 


i^-yjT'C^ + yjf'dyj, 


(71) 


13 






where o.j,(3j > — 1 with the uniform measure as the special case aj = fij = 0 , the corresponding 
orthonormal polynomials are the Jacobi polynomials given by the Rodrigues’ formula 




'’-w = n + nr”’ (fej - - viT-a +•/j)*), (72) 


i>i 


where 


0!j,l3j _ 

Cui ■ — 


(2i7j + cKj + fjj + 1 ) i/j! r(i7j + cxj + fdj + 1 ) r(Q(j + 1 ) r(/3j + 1 ) 


r(i7jr + Q(j +1) r(i/j -|- Pj -|- 1 ) r(Q;j + fdj + 2 ) 


(73) 




with the convention that Cq^’ ” = 1 for any aj,l3j > —1. Using integration by parts analogously to 

(1 - y^P 


(63); for the Jacobi coefficients Uy we then obtain 
Uy=[ u{y) Py{y) dfi{y) = Oy [ ^d''u{y) 


lu 


'u 


i>i 


2'^i 


-dpy), 


ay Cy, 

i>i 


aj ,/5j 


(74) 


For p = {pj)j>i satisfying the assumptions of Theorem 3.1 proceeding exactly as in (64) we arrive 
at 

'^{dpp''\\uy\\v)‘^ < oo. (75) 




Provided that the sequences {aj)j>i, iP )j>i are such that Cy-’ ” < Qpj) for some polynomial Q 
independent of j (which holds, for instance, when both sequences are constant), the statement of 
Corollary 3.1 thus holds also for ( 11711/11 v')yg_;r- 

Remark 3.4 The same remarks as given at the end of §2 apply to the generalization of Theorem 


3.1\ Corollary \3.1\ and Remark 3.3 to other types of linear PDFs with affine dependence in the 
coefficients. 


4 Examples 

In this section, we compare the summability properties obtained with the approach in the present 
paper to those obtained with the analysis in for various types of (V'jOjM- We show that this 
approach gives an improvement on Theorem El depending on the particular structure of the 
supports of the i^j. 

4.1 Finitely overlapping supports 

We first consider families Pj)j>i of functions with hnitely overlapping supports, that is, such that 
for any x € D, 

#{j : V',(^)^0}<M, (76) 

for some fixed M > 0, where ffS denotes the cardinality of the set S. The case M = 1 corresponds 
to disjoint supports, such as the family of characteristic functions 'ijjj = bjXoj with some normalizing 
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factor bj, when {Dj)j>i is a partition of D. Another example with M > 1 is the set of Lagrange 
hnite element basis functions of a given order A: > 1, associated to a conforming simplicial partition 
of D. 


Assuming (UEA*), we then dehne a weight sequence {pj)j>i by 

hmin(l 


Pi 1 + 


2M||'(/;j: ll^oo 


(77) 


With such a choice, we hnd that (39) holds since, for all x G D, 


^^ < 5a{x), (78) 


i>i i>i 

with 5 := <1. As a consequence of Corollaries 


^ ircjiiL” 

and 


2.1 


3.1 


we obtain the following result. 


Corollary 4.1 Assume that is a family of functions with finitely overlapping supports, and 

that I^UEA *) holds. If (||V’illL°°)i>i ^ ^"^(N) for some 0 < q < oo and if 0 < p < 2 is such that 


q = qip) := then (||ti/||y)i/ej' and (||u,y||y)i/gj- belong to ^{IF). In particular, best n-term 


2-p 


Legendre approximations converge in L‘^{U, V, p) with rate n * where s = - — 4 = -. 


As already observed we always have q{p) > p which shows that there is in this case a signihcant 
improvement between the summability properties of {\\'4’j\\L°°)j>i and those of (||t,y||y)j^gj- and 
(||tfi/||\/)i/ej^, in contrast to Theorem 0 Note in particular that the latter would lead to the 
weaker conclusion that best n-term Legendre approximations converge in iP'fU, V, p) with rate n“® 
where s = - — i instead of -. 

q z q 

We next give a specihc example which shows that, for the Taylor coefficients, this new result is 
in fact sharp. In this example, we let D :=]0,1[ and a = 1, and we consider a sequence of 

disjoint intervals Dj =\lj,rj[(Z D. Let mj := ^{Ij + rj) be the midpoint of Dj and ijjj := bjX[i.^rnj] 
with {bj)j>i G .^'^(N), where q = q{p) := for some 0 < p < 2. We denote by 

hj{x) := maxjO, 1 — 2|x — mj\/\Dj\^, (79) 


the hat function on Dj centered at ruj with hj{mj) = 1. We hx a sequence {cj)j>i such that 
^i>i ^ ^ choose the right hand side / = —(X]j>i G V so that 


to-Y 


Cj hj. 


i>i 


(80) 


The condition on the Cj ensures that tg £ Foi' the particular u = ej, the Taylor coefficients 
satisfy 


D 


t^.v dx = — ipjtoV dx , V ^ V . 

^ Jd 

Testing this with v = hj, by the Cauchy-Schwarz inequality we obtain 


Iti.P dx 


ID 


D 


D in.12 

up 


> \Dj\bjC_ 


I n 12 ’ 

Uj\ 


(81) 


(82) 
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and hence 


(83) 


ll^e 


If > 


bjCj 

V\Dj 


In view of the requirements {bj )j>i £ £'?(N) and {cj/^/\Dj\)£ £^(N), these sequences can be 
chosen to ensure (||tej || V')j>i ^ ^^(N) for any p < p. 


4.2 Arbitrary supports 


In the general case where the supports of the '0^ are arbitrary, in particular for globally supported 
functions, the approach based on Theorem 2.1 does not bring any specihc improvement over Theorem 
o (which can be derived from it, as observed in Remark |2.1[ ). One way to see this is to observe 
that in certain situations, the latter can already be sharp. 

Consider for example the case of constant 'ijjj = bj with {bj)j>i £ ^^(N) and {hj)j>i ^ .^^(N) for 
any p < p, and a = 1. For such tpj, one has 


u{y) = 


1 


1 + Ej>i Vjbj 
so that the Taylor coefficients are explicitly given by 


?x( 0 ). 


(84) 


A = (-l)l"l^6%. (85) 

In particular, one has UtejUv' = which shows that (||ti/||y)i/ej^ ^ for any p < p. A 

similar, yet more technical, computation shows that the same holds for the Legendre coefficients. 
Therefore, in this case of completely overlapping supports of the the proposed new bounds 
cannot give an improvement over Theorem o 


4.3 Wavelets 

Let us now turn to the case of diffusion coefficients parametrized by a wavelet basis, that is, 

a{y) = a + '^yx'ilJx, (86) 

A 


where 

IIV'aIIl-= Q := C2-“', |A|=/, (87) 

for some a > 0, as discussed in the introduction. Note that, when ordering the wavelet basis from 
coarse to hne scales, the resulting system (V’j)j>i has then the algebraic behaviour 

llV’ilU” ~ (88) 


Assuming (UFA*), and for an arbitrary 0 < /3 < a we dehne a weight sequence {px) by 


Pa := 1 + 


Omin(l ^) 
2C'M(1 - 2/5-^ 




(89) 
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where M and C are the constants in (14) and (87). With such a choice, we find that (39) holds 
since, for all x G H, 

< flaw + E ^ ~ 


2CM{1 - 2/5-“) 


< Sa{x), (90) 


where 6 := < 1. After the same reordering as for the wavelet basis, we find that 


Pj 


nP/d 


(91) 


Therefore, as a consequence of Corollaries 2.1 and 3.1 we obtain the following result 


Corollary 4.2 Assume that {'4’j)j>i is a wavelet basis with normalization (81) and that (UEA*) 
holds. If i\\'4’j\\L’^)j>i £ f"^(N) for some q < q{p) := then i\\tu\\v)u£T {\\uu\\v)u£T belong 
to (P{F). In particular, best n-term Legendre approximations converge in Lf‘{U, V, p) with rate n”® 
for all s < -. 

J q 


As already mentioned, if we use sufficiently smooth wavelets, the decay property (87) is equiva¬ 


lent to the property that a{y) is in the Besov space (D)), which for non-integer a coincides 


with the Holder space C°‘{D), for all y £ U. Thus, we also infer from Corollary 3.1 that if this holds 


for some a > 0, best n-term Legendre approximations converge in Lf{U, V, p) with rate n ^ for all 
s < ajd. 

5 Numerical illustrations 


In the following numerical tests, we consider three different cases of parametrized diffusion problems 
on D =]0,1[. In each of these cases, the parameter dependence is expressed in terms of a different 
representative type of function system i'4’j)j>i that corresponds to one of the three types of examples 
considered in Section HI 

We give results both for Taylor and for Legendre coefficients. The Taylor coefficients are com¬ 
puted using the alternating greedy Taylor algorithm described in [TJ §7.2]. For the Legendre coeffi¬ 
cients we use an adaptive stochastic Galerkin method using the approach described in |12| . where 
multi-index sets are refined based on approximate evaluation of the residual of the continuous prob¬ 
lem. Both methods are guaranteed to converge, but in the Taylor case one needs to resort to a 
heuristic stopping criterion. For each of the tests, a suitably adjusted single fixed finite element 
discretization using P2 elements is used for the spatial dependence. 

In all of the following examples, we take d := 1. We also fix 9 g]0, 1[ and choose the functions 


ifj in the examples such that this 6 satisfies (29). 


In the following discussion of the methodology for evaluating the results, we always refer only to 
the Taylor coefficients (t,y)jygjr, with the understanding that the same considerations apply to the 
Legendre coefficients {ui,)u£T- In each of the tests, the aim is to numerically estimate the limiting 
summability exponent 

p:=inf{p>0: {\\ty\\v)u(iT ^ P{IF)] ■ (92) 
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To this end, we introduce the decreasing rearrangement (tjj)n>i of (||tj^||y),ygj-. Then G 

P’{F) implies that for some C > 0, one has t* < and conversely, if t* < for some 

C,q > 0, then (||t,y||y)i,gj- G l’P{F) for any p > q. As a consequence, 

1 

p = inf|p>0: supnpt*<oo|, (93) 

nSN 

or in other words, p can be determined from the asymptotic decay rate of the values t*. As estimates 
for the largest s > 0 such that sup„ is finite, we consider the values 


Si := log2(t2*-i) - 


(94) 


for i = 1,2,.... In view of (93), if the sequence (t)j)n>i decays asymptotically at an algebraic rate, 
for sufficiently large i we can thus expect Si to approximate p~^. 


5.1 Parametrization by disjoint inclusions 

In the first test, we choose a family {Dj}j>i of disjoint open intervals in D and a /3 > 0, and define 

^.:=0j-PXDy (95) 


Note that although this does not enter into any of the decay estimates available at this point, the 
concrete example in |4.1| suggests that the decay of the inclusion sizes \Dj\ has an impact on the 
summability of (||f,y||y)[/gj- and (||tt;y||y)j,gj-. Indeed, it can also observed numerically that faster 
decay of \Dj \ leads to improved summability. To remove this effect in the tests, we therefore choose 
Dj such that the decay of \Dj\ is as slow as possible while still allowing Dj to partition D. To this 
end, we define 

3 

xo := 0 Xj := log“^(l + fc), j > 1, (96) 

k=l 


with c such that limj_>.oo Xj = 1, and set Dj =\xj-i,Xj[. Since 
by Corollary 


4.1 


'j\\L°°)j>i £ ^"^(N) for all q> 


we expect that (||t;y||y)!yeJ' (ll^!^l|y)i.'eJ' belong to IIP{F) for any p > (/3+ 4) 


u-i 


The values of the decreasing rearrangements of these sequences for /3 = ij,l,2, where in each 
case 9 = ^, are compared in Figure In Table the empirically determined decay rates are 
compared to the theoretical prediction for p~^. We observe almost the same decay behavior for 
Taylor and Legendre coefficients, and in each case the empirical rates indeed approach p~^. 


5.2 Parametrization by a Fourier expansion 

We next consider a parametrization by the globally supported Fourier basis 

V’j(x) := sin(j7 rx), 


(97) 


for some /3 > 1, with the normalization constant c := ^) ^. We thus have 

due to Theorem 


£P(N) for all p > In view of the discussion in ^4.2 
and (||n[y||y)jygj- belong to (P{F) for such p. 


1.1 


we expect that (||t,y||y)i,g_ 7 r 
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Figure 1: Disjoint inclusions: ordered norms of Taylor coefficients ti, and Legendre coefficients u,^, 
for 0 = ^ and /3 = ^, 1, 2. 




Taylor 


Legendre 


II 

to 

/3 = 1 

/?=5 

II 

to 

/? = ! 

/3 = 5 

S6 

2.563 

1.730 

1.225 

2.476 

1.789 

1.302 

S7 

2.708 

1.731 

1.274 

2.578 

1.786 

1.235 

S8 

2.481 

1.726 

1.211 

2.601 

1.701 

1.212 

sg 

2.574 

1.706 

1.235 

2.514 

1.661 

1.200 

Sio 

2.439 

1.650 

1.196 

2.543 

1.660 

1.169 

5ii 

2.477 

1.643 

1.175 

2.507 

1.642 

1.160 

p-^ 

2.500 

1.500 

1.000 

2.500 

1.500 

1.000 


Table 1: Disjoint inclusions: decay rates of coefficient norms, with s* as in (94), compared to limiting 
value p~^ = /3 + I expected by Corollary 


4.1 




Figure 2: Fourier expansion: ordered norms of Taylor coefficients and Legendre coefficients Uiy, 
for 0 = I and /3 = 1.25,1.5, 2. 
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Taylor 



Legendre 



(3 = 2 

II 

13 = 1.25 

II 

II 

CJl 

= 1.25 

^6 

1.452 

1.165 

1.250 

1.593 

1.294 

1.250 

S7 

1.619 

1.320 

1.092 

1.682 

1.353 

1.154 

ss 

1.495 

1.278 

1.147 

1.597 

1.337 

1.192 

sg 

1.515 

1.257 

1.141 

1.632 

1.338 

1.187 

SlO 

1.533 

1.270 

1.143 

1.637 

1.341 

1.173 

Sii 

1.515 

1.258 

1.143 

1.639 

1.327 

1.191 

p-^ 

2.000 

1.500 

1.250 

2.000 

1.500 

1.250 


Table 2: Fourier expansion: decay rates of coefficient norms, with Si as in (94), compared to limiting 
value p~^ = j3 expected by Theorem 0 




II 

to 

1 

e = 2-3 

e = 2-3 

•^6 

1.593 

1.876 

2.000 

S7 

1.682 

1.767 

1.872 

S8 

1.597 

1.822 

1.908 

sg 

1.632 

1.813 

1.905 

SlO 

1.637 

1.813 

1.898 

5ii 

1.639 

1.814 

1.921 


2.000 

2.000 

2.000 


Figure 3: Fourier expansion: ordered norms of Legendre coefficients for /3 = 2 and 0 = 
and corresponding decay rates. 


The results for 0 ^ and j3 = 1.25,1.5,2 are shown in Figure and Table Here we observe 

that especially for larger values of /3, the empirically observed rates Si do not come very close to 
the theoretically guaranteed limiting value p~^ within the considered range of coefficients. This 
indicates that the asymptotic behavior emerges only very late in the expansions. 

Note that this observation is consistent with the results obtained in [12], where the above 
example (97) with (3 = 2 and 0 = ^ is considered as a numerical test. There, a decay rate close to 
1 is observed for the error of a Legendre expansion (with fixed spatial grid), corresponding to a 
decay rate of the coefficient norms close to 1.5 as obtained here. 

It turns out that the observed decay rates are in fact also influenced by the value of 9: as 
shown in Figure for smaller values of 6, the Si are closer to the limiting value already within the 
considered range. 
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Figure 4: Haar expansion: ordered norms of Taylor coefficients ti, and Legendre coefficients Ui,, for 
0 = ^ and a = ^, 1, 2. 


5.3 Parametrization by a Haar wavelet expansion 


As a final example, we return to the wavelet parametrization of a with a levelwise decay (87). Here 
we use the Haar wavelet, generated from h:=XjQij — such that 


'ipx{x) = cih{2^x — k), X = {l,k), / > 0, k = 0,. 


. 2 ' - 1 , 


and we set 


Q := 0(1-2-“) 2 


—al 


(98) 


(99) 


for a hxed a > 0. Since, after reordering, we have 
for all g > ^, by Corollary 


.jil^oo ~ j ^ and therefore (ll^illL-).>ie^nN) 
we expect that (||ti/||v')i./ej' and (||u,y||y),ygj- belong to (P{F) for any 


C-i 


4.2 


p > (a+ i)- 

The results for 0 = ^ and a = 1, 2 are given in Figureand Table We again observe very 

similar decay for Taylor and Legendre expansions. Similarly to the observations made in Section 
5.2 for larger a, the empirical rates Si do not come very close to the expected asymptotic limit p~^ 


within the considered range of coefficients. As shown in Figure the decay rates again approach 


p ^ more quickly for smaller values of 0. 


In summary, these numerical results support the conjecture that the summability estimate in 
Corollary |4 .2 1 for wavelet expansions is in fact sharp, similarly to Corollary [4T] for disjoint inclusions, 
whose sharpness we have established by an analytical example. 
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Taylor 


Legendre 


0 = 2 

a = 1 

0=1 

0 = 2 

a = 1 

0=1 

^6 

1.450 

1.301 

0.927 

1.853 

1.165 

0.947 

S7 
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II 

to 

1 
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S6 
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S7 
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